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ABSTRACT 


In this paper a unified presentation is made on the results of var- 
ious investigators on the properties of random communication networks. 
These results are interpreted in such a way that the properties may be 
determined by using a digital computer with the application of the 
Monte Carlo method. The computer program is written and tested. Re- 


sults for some networks are compared with theoretical values. 
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I. RANDOM COMMUNICATION NETWORK 


In this paper the results of various investigators on the properties 
of random communication networks are reviewed. From these results a 
mathematical model is represented for a general random communication 


network. 


A. THE CONCEPT OF A RANDOM COMMUNICATION NETWORK 

A communication network is an aggregate of message centers that at- 
tempt to transfer information to one another over a wide variety of 
channels. Themessage centers could be fixed as the headquarters of a 
regional military command, mobile as ships, aircrafts, or satellites. 
The channels could be a radio link, a telephone line, or highway. In 
normal conditions, the topology of the network is known with certainty: 
the positions of ships at sea, the links between the headquarters with 
different units. Each center knows its relative position with respect 
to the other centers exactly. If a message is sent from a central 
command, the authorities know exactly how many and which units will re- 
ceive it. The network, in this case, is called a deterministic 
communication network. But during the hostilities, some units could be 
— by the enemy without even having the chance to send the last 
word, some radio link could be jammed by the enemy, some telephone lines 
or some roads could be damaged without having any means to determine the 
extent of the destruction. Back at the headquarters, the authorities 
are faced with the uncertainty of the existence stage different units 
and channels, the effectiveness of an order sent to the battlefield is 
not as high as before. It is only probable that such unit gets the 
message. The communication network becomes a highly probabilistic, 


random netwrok in the event of a war. 
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One finds the same kind of randomness in many other branches of 
research; for example, the contagion of a disease among a population, 
an hereditary defect in a species, the spread of a stimulus in the 
nervous system, or the spread of information in society, to name a few. 
All these phenomena are very similar to the random communication net- 
work in that they have: ; 

a. an originator which originates a specific message 

b. several recipients which convey the received 'message' to the 
others via a medium (atmosphere, voice, physical contact, etc.) 

c. the uncertainty of contact between two individuals. 

All these phenomena could be then studied under the same topic of 
the random communication network. Researches in these phenomena give 
much insight in the problem of propagation message through a random 
network. 

As it is said, the communication between = message centers A and 
B could be assumed only with a certain probability. The latter measure 
is a function of many parameters which are for example, the time t 
when the message is sent, the time T since the message is received at 
ón A, the distance d between stations, the electronic equipment of 


each station, to cite a few. Therefore one could write, in general: 
Puobe(econcact A to by) = pe(A, B, t,T,d,.. .). 


Due to the complexity of the problem, without reasonable assumptions 
about the parameters to simplify the problem, one could not solve it. 
Following is the brief survey of the contributions of various workers. 


Each of them has solved partially the problem with different assumptions 





about the parameters which determine the probability of contact between 


two stations, two individuals, whatever is the case. 


B. THEORETICAL DEVELOPMENT 

Bailey [1] studied the spread of an epidemic in an homogeneous 
population with size n where the probability of one new infection tak- 
ing place in the interval dt is given by 

Prob (one infection in dt) = y (n-y+1)dt, 

where y = the number of susceptibles. 
From the above equation, the mean m(t) of the number of infected people 
and the epidemic curve, (i.e., the rate of change with respect to time 


of the mean of the number of infected people) dm(t) are obtained to be 
dt 


m(t) = e"20tc810t - 234) + e718%4410t - 902) + . 


dm(t) 


= е-19% 81006 - 3156) « e-18t(79,380t - 20650) + . 
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Landau and Rapoport [2] assumed the probability of conta ct between 
any two individuals is the same for any pair and the probability of 
transmission of the disease depends both on the time t of the whole 
process and also ont, the time since the particular affected indivi- 
dual acquired the disease. Thus, 

Prob (one infection) = p(t,T) 

The differential equations giving the number of individuals who become 


affected and the rate are: 
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where N - number of individuals in the population 
x(t) = total of individuals who have become affected up to time t 
хо = number of individuals who become affected at t = 0, the 
initial time of the process 
A =frequency of contacts, which is assumed constant. . More pre- 
Cisely dx(t)[N - x(t)] is the number of contacts per unit time 


between affected and unaffected individuals. 


This model is more realistic than that of Bailey because it is true 
that in the case of a disease epidemic, the infectiousness of the path- 
ogenic virus may increase or decrease with the time of the process and 
the time since it infects the body of the individual. 

The equations (1,1) and (1.2) are difficult to solve in the general 
case, therefore the investigators made some assumptions and solve them 
for the following — 

а. The probability p(t,t) is a function of t alone 

b. The probability p(t,t) едиа15 е-КС 

c. The probability p(t,t) is function of t only 

d. The probability p(t) is constant in a finite interval 

e. There is a lag in the transmission of contagion. 

Landahl [3] treated the spread of some information as a flow of 
'particles' which execute random motions over a population of individuals 


and which may multiply or disappear. Equations are derived for the 
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density of these particles, and the number of individuals through 
which the 'particles' have passed, is calculated. The results are ap- 
plied to the following cases: 

a. Uniform spatial distribution with multiplication factor decrea- 
Sing with time because of loss of interest. 

b. Multiplication factor is constant, but the rate of spread de- 
creases with multiple hearings. 

c. One dimensional region with a small starting region with or 
without an absorbing barrier. An absorbing barrier corresponds to the 


case where the individuals listen to the information but do not repeat 


it. 

d. Two-dimensional region with absorbing barrier. 

e. Continuous sources of information within a small region in one 
dimension. 


f. Uniform spatial distribution in which individuals do not 
respond to more than one hearing. 

Prihar [4] studied a mobile communication network where the trans- 
mitting and receiving stations are in motion. He used the classical 
Kinetic Theory of gases to compute the expected total number of vehi- 
Bios! comta per unit time in two cases: 

a. The search for contact is continuous and the antenna is omni- 
directional. 

b. The search is intermittent and the antenna is directional. 

Finally, Mattei [5] assumed the probability of contact between any 
pair of stations is constant throughout the network. 

Prob (contact between any pair) = p = constant. 


Several formulas were derived in closed forms, 
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The mean number of contacted nodes at the lst step is 
2(1) = (n - р. 

The mean number of contacted nodes at the 2nd step is 
102) = (п - 1) (1 - (1 - р2п-21. 


The marginal probability of the number of new nodes at each step is 


y (0) 1 k-2) 


Prob[Z(k) = z(k)] = Y E 
z(1) - 0 z(2)-0 2 (КӘ) = 


(96) EOS [a2] Q) 


2(2) 2 УС) 


] 


Ua) fil = са [а 


2 (2) 


со [1 - q2(k-1)]2(k) EMEN 


SK 
y(k) =n - z (i) q = 1-p (1.5) 
1=0 


where n = number of nodes in the network 


k = order of step 


C. THE PROPOSED MATHEMATICAL MODEL 

Although with the model of Mattei, one could derive many quantities 
of the network in closed and neat forms, the model is not realistic in 
real life as pointed out in the beginning of this paper. The probabil- 


ity of contact between any pair of stations is far from constant. It 
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depends on their distances to the center of attack, the size of the 
electronic equipments, the geographical positions. In the event of a 
war, the capital, the missile silos, the industrial centers are more 
vulnerable than a submarine under the sea or a communication satellite 
In orbit. 

In this paper, a model, in which the probability of contact between 
any pair of stations could be varied at will, could be assigned to each 
link, once all circumstances have been taken into account, is proposed 
and solved. The solution is obtained by simulation on the digital 
computer. The computer program allows the user to assign the probabil- 


ity of contact between all pairs of stations. 


ie 





ІІ. SIMULATION ON A DIGITAL COMPUTER 


The aim of this chapter is to obtain, by simulation on 
the digital computer using the Monte Carlo method, the de- 
sired quantities of a general random communication network, 


in which the link probabilities could be all distinct. 


А. THE MONTE CARLO METHOD 

The Monte Carlo method consists of solving various 
problems by means of the construction of some random process 
ror each such problem, with the parameters of the protess 
equal to the required quantities of the problem. These 
quantities are then determined by means of the observations 
of the random process and the computation of its statistiCal 
characteristics which are approximately equal to the re- 
quired parameters. 

For example, the required quantity x might be the math- 
ematical expectation EN) of a certain random . The Monte 
Carlo method for determining the approximate value of the 
quantity x consists of an N-fold sampling of the vane of 
the variable in a MGE: of independent tests, (that is, 


1,12, ши? ) and the computation of their mean value 


И = DAS ET ы ыз в Ум 
М | 


12 







— адан i ч EE -- 
a, p ML, Lu : 


— > 
Å 





o 


ncm ———— |. н D 
е. — an с a m 
a mom m c ж | 44 s ome 
då >» e wq o c——p» «к © O 

ж å bly c— o Â ID 





If the number of tests is large, one has 


It is appropriate to ask the question why one does not 
try to derive an analytical expression for the quantity x 
from the properties of the problem instead of designing 
a random process which imitates the reality. The answer 
is that many problems in the physical world are so complex 
that they do not have closed formula for the desired quan- 
tities. 

Suppose, for example, that the'required quantity x 
is the definite integral of a function f(y) taken over 


the interval (a, b), that is 





Fig. 2.1. Illustrating equation (2.2). 
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If the expression of f(y) is not simple (or even worse, 
if the function f(y) does not have a closed formula but 
is known only point by point), x could not be evaluated 
analytically. The Monte Carlo methOd may be used as will 
be described next. 

Suppose we have a device called a random number gen- 
erator which generates two random variables E and | 
uniformly distributed cver the intervals lag bp)” and Vorme), 


respectively. The Compound probability density function 


is given by 


1 
с(Ь-а) if açęyç¿b and 0O<z<c 
Pg sh (y,2) = 
0 otherwise 


The pair (En) will be generated N times and each time, the 


condition 


zs De е (29 


will be tested. If condition (2.1) holds, the point 
(Е, Ni) in Fig. 2.1 is inside or on the boundary of the 
cOmputed area. Let N be the number of times when this 


condition holds in N tests, we have 


b . / 
; _ lim (N (2.2) 
х = f(y)dy = c(b-a) (2) 
J N-»o^ М 
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In order to obtain good results, the number of tests 
must be large, say, in the range between 1,000 and 100,000. 
This is only possible with the aid of digital computers. 

The problem of a random communication network where 
the probability of survival of each link is assumed constant 
has been solved adequately by Mattei [5]. If one changes 
this assumption and lets the said probability vary to fit 
certain Conditions in reality, the problem becomes Umtrac- 
table, then the only feasible solution is to use a 
numerical method such as the Monte Carlo technique. The 
required quantities are: the average number of contracted 
nodes, the average number of newly contacted nodes at each 
step, the terminal reliability, the probability distribution 
of new nodes at each step. The method for obtaining these 


quantities will be described in the following section. 


В. THE FLOWCHART OF THE SIMULATION PROCESS 

Suppose that one has a net—work consisted of n stations 
connected by a certain number of links, it is a determinis- 
tic network in the sense that, if a message is sent by a 
ЕЕ са to adjacent stations which relay it to the others, 
one knows exactly how many and what stations have received 
it. Suppose a disaster, which could be a nuclear attack, 
am earthquake, etc., happens to the network. Before any 
investigation of the damage could be alos one does not 
know which links have survived. If a message is sent, it 
is difficult to acertain how many and which stations will 


gct the message. The best one сома do ise to assign to 
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each link a certain probability of survival. Let DER be 
the probability of communication between station i and sta- 
tion j, the probabilities between all pairs of stations are 


best described by an n x n matrix as 


E RENE ші 
ie hd ee IR 
РЕ 
Р. 
Poy Py -- -Pn, n-1 0 


= 


where Р 15 called the matrix of probabilities. 


Since one 1S not interested in the case where a 


station sends back to iteself a received 
message, all diagonal elements of P are zero. From a set 
n(n-1) networks. Бъри 2. 2 


of n stations, one could form 2 


shows two of them for n = 5. 


e M ж 


=, 


© © 


Fig. 2.2. Samples of random network 
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If the number of stations increases, say n = 21 stations, 
the population of possible networks is 2420 ха E 29x 
6°. which is large. Suppose that one investigates the 
propagation of message via the originator which is the 
station 1. Starting with this station, one generates a 


random variable r uniformly distributed between O and l, 


and check the condition 


A 


O 9 
— ! 





© 
o | 
{ 
2 | 
i 
| - | 
| | 
| | 
| @. | 
| I 
| | ЄР | | 
| | (4) 
step O зрер I Step 2 Stcip э 


2(0)-1 2(1)=2 2(2)=4 2(3)-2 2(4)-0 


Fig. 2.3. A tree of propagation 
PAR В ( 25989) 


If this condition holds, there is a contact between station 
] and station 2 and one increases the random variable z(1), 


which represents the number of newly contacted nodes at 
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ети, бу боле шла Е. One generates again a random number г, 


fests comdltiven (2.5) for P and 3mcreases z(1) by one 


1,3 
unit if condition (2.3) holds, and so on until one arrives 
at the last station, which is station 21. One repeats the 
same procedure with each station of z(1) relative to the 
remaining network, the number of stations contacted at this 
second step is designated z(2). One repeats the same pro- 
cedure with each station of 2(2). The procedure stops at 
the step when one could not contact any more stations. 

As an example, consider the network in Fig 2.3, stations 
5 and 17 have been contacted at step 1 so that z(1) = 2. 

The stations which have received the message are discarded 
from the network, these are stations 1, 5, 17. Next, to 

each station contacted at step 1, the same procedure made 

for station 1, is repeated. In the example, stations 2, 6, 

8 and station 14 have been contacted respectively by station 
5 апа station 17 at step 2, z(2) = 4. Stations 4 and 13 

have been contacted by station 6 at step 3, z(3) = 2. The 
propagation of the message is over at step 4 because stations 
4 and 13 cannot carry the message any further. 

The tree in Fig. 2.3 summarizes the process in various 
steps. At any time the computer needs to keep the infor- 
mation of only one tree, which saves a lot of memory 
Storage. To start the next sampling, one just has to erase 
all information of the first sampling run except, of course, 
Ее resmMts z(1), 2(2), . . . . The general procedure for 


the second sampling is the same as for the first. 
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After K samplings, the average number of new nodes at 


each step is given by 


ON ОКЕ NNN 2 С) 
_ ЖЕЛІ ae NN + 2 12) 
ee > a 

21 (КМАХ) + 2. (КМАХ) т + 2, (КМАХ) 


where KMAX = the highest order of steps in the tree. 


The average number of contacted nodes is 


ENN = о. 


Се 2.2 + 2. (КМАХ) 
A о + Z, (KMAX) 
ШӨ = 2 (1) + 2, (2) +... + 2, (КМАХ) 
Х = = Е =. ВЕ 
К 
= X(1) + X(2) + . . . + X(KMAX) 
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where X(i) is the average number of contacted nodes up to 
sep i. 


The Kth terminal reliability is commonly defined as 


Z(K) 


r(K) = 1 





3 


and the weak connectivity is defined as 


Х(КМАХ) 


ү > == 


Definition. In this paper the link density of the network 
is defined as the ratio of the number of links between sta- 
tions to the number of emitting stations. 

During the process of simulation, the probability dis- 
tribution of new nodes at each step is also calculated and 
displayed in an n x n matrix whose rows represent the num- 
ber of steps and columns represent the number of new nodes 


as follows: 


# of contacted stations 


0 1 2 3 å f ы n 


order of step 


n 


Fig. 2.4. Probability distribution of the number of con- 
Bacted stations at each step. 
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A flowchart of the simulation is drawn in detail in 
fee. 255. The list of the variables used in the computer 


program is described as follows. 


IRMAIN = vector whose components have the value 1 or 0. 
If the ith component is 1, it means station i 
has yet been contacted, 0 means the station 
ith has been contacted. 

к 

NEW = vector whose components have values 1 or 0. If 
the ith component is 1, it means the station 1 
has just been contacted during the preceeding 


step. 


TROPY = vector with 1 ог аз сопропеде ад the ith 
component is 1, it means the station i has been 
contacted during the current step. Before to 
go to the next step, one has to make; NEW = 


ICOPY. 


ZVEC = vector whose ith component represents the num- 
ber of newly contacted stations at step (i - 1). 
XVEC = vector whose ith component represents the aver- 
age number of contacted nodes at step (1 - 1). 
GVEC = vector whose ith component represents the 


(i - l)th terminal reliability. 
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RELIAB 


RNGE 


PRDIST 


PROB 


RLIAB 


AVEC 


vector whose ith component represents the. 


(i - 1) terminal reliability. 


vector whose ith component represents the prob- 


ability distribution of new nodes. 


vector whose components determine the x, y 
ranges of the plot in the subroutine UTPLOT to 


plot the average number of new nodes. 


Square matrix whose (i, j) component represents 
the probability that the number of contacted 


nodes iS j at step (i - l). 
square matrix whose (i, j) component represents 
the probability of communication from station 


i to station j. - 


Square matrix whose ith row represents the ter- 


minal reliability for the ith case (different 


probabilities, different standard deviations, 


ECC nae 


rectangular matrix whose ith row represents the 
average message propagation for the ith case 
(different probabilities, different standard 


deviations, etc...). 
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KSTEP 


KHECK 


КНЕСКІ 


5 МКО 


ZNODES 


KMAX 


KKMAX 


CONECT 


NDEMIT 


LABEL 


MC 


ТХ 


INC PersUhNcehwrepresents the order of tife step. 
integer whose value is 1 or 0 used to initialize 


ICOPY to (0) or to compute the number of emitting 


stations. 
number of contacted nodes up to step K. 
number of contacted nodes during step K. 


highest order of step attained during the pro- 


cess of simulation. 


KSTEP + 2 used in the computation of PRDIST 


minor variable. 


represents the weak connectivity of the network. 


the number of emitting stations during one sim- 


ulation. 


used in the subroutine DRAW. 


integer used in the subroutine DRAW. 


integer used in the subroutine OVFLOW and RANDOM. 
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Капа Z 


Р апа РХ 


the average number of contacted stations at 
the lst step in the model of Mattei or defined 


as n/A, standard deviation in the general model. 


the standard deviation of destruction in the 


general model, defined as: S = n/A. 
Normalized variables, defined as X = I/S and 
ель. 


the number of tests in the process of simulation. 


outputs of the subroutine NDTR which are areas 


under the normal curve. 


The following subroutines of the NPGS computer center library 


are used in the program. 


UTPLOT 


DRAW 


OVFLOW, 


to plot the average number of new stations curve. 


to plot the probability distribution of con- 


nected stations, the terminal reliability and 


the average number of contacted stations. 


RANDOM = used to get a random variable (0,1). 
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NDTR | = to obtain the area under the normal curve as 
probability of communication between two sta- 


tions. 


To investigate the general model where the destruction 


is assumed normally distributed around the center of attack, 


block(2) has to be replaced by block(l). 
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SIMULATION OF A RANDOM COMMUNICATION NETWORK 














ZVEC 





IRMAIN - (11. . . 1) (10 0) 


NEW = (00... 0) A *. . 0) 
ICOPY = (00 . . . 0)  GVEC= (0... 0) 
KSTEP = 0; KHECK = 0, КНЕСК1`= 0 
SUMND = 1.0, ZNODES = 0.0, KMAX = 0 
KKMAX = 0, CONECT = 0.0, DNSITY = 0.0 
NDEMIT = 0, МС = 1, IX = 102530469 


RANGE» = (20.0, 0.0, 2140, 040), К = 38,000 


| CALL OVFLOW | 


pM ыы | 






КЪТАВ <- (0) 









СОС -MOM 2 v l N) 


XXVEC < (0) 


Fig. 2.5. The Flowchart of the simulation. 
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STOP |-- 
p. | PROPAGATION  -— 
E 





PLOT: TERMINAL | 


RELIABILITY AND, 


CURVES p 


BLOCK (1) 
Dec S UO PRO (1J) = 
S «- 21/A PROB (J,I) 
X « I/S PROB (I,J) = 0 
САШ, МОТА (I= J) 
2 < 2/5 
CALL NDTR 


PROB (1,J) < (PX - 0.5)MP- 0.5) 


-— M 00 A 





ЕТР A 


LG 


“ОО 






CURVE 


С vem 


BLOCK (2) 


— i 
| 


А = А + 1.0 


А 
PROB © (A 
S 


PRDIST < (0) 


Lo e e a dil e MR E A 







PRDIST (1,2) = M 
IRMAIN (1) = 0 
NEW (1) = 1 


ICOPY (1) = 1 


continued. 
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DNSITY < 0.0 
CONECT < 0.0 
ZVEC < (0) 


RELIAB < (0) 


reu» 
mu 


DNSITY 
MESSAGE, SPEED 
PROPAGATION 


PRDIST 


ZVEC A 


MUN | 
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| 
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CALCULATION OF: 
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PROPAGATION SPEED, 
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С. COMPUTER PROGRAM 

The computer program is written to analyze a network of 21 sta- 
tions, however this number could be increased at will. First one has 
initialize the matrix of probabilities of communication between sta- 
tions to desired values. The probability distribution function could 
be discrete or continuous. 

For example, in the case where the probability of communication 


between any pair of stations is constant, it suffices to write 
Table 2.1. Fortran statements for the model where p is constant 


DO 16 I-1,N 
| DO 15 J=1,N | 
IF (I. EQ. J) GO TO 14 | 
PROB (I,J) = A/(N - 1) 
GO TO 15 
14 PROB (I,J) = 0.0 
15 CONTINUE 


16 CONTINUE 


EE Na 


If one wishes to analyze a network where the destruction is normally 
distributed around the center of a nuclear attack, it suffices to 


replace table 2.1 by the following table 
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Table 2.2. Fortran statement for a normal model 
$ = 21/А 
DO 16 I = L,N 
X = I/S 
CALL NDTR (X, PX, D) 
DO 15 J = I, N 
FF (I. EQF J) GO TO 14 
2 = J/S 


CALL NDTR (Z, P, D) 


РОВ (I,J) = (PX - 0.5) . (P - 0.5) 
PROB (J,1I) = PROB (I,J) 
14 PROB (I,J) = 0.0 


15 CONTINUE 

16 CONTINUE 
The computer program could be used to analyze networks of any size 
N by redimensioning the variable names which are vector arrays, matrix 
arrays to the size of the network and enter the data card according 
to ee changes. The program prints out the link density and the weak 
connectivity of the network, as well as the curves of the average 
number of contacted stations, the probability distribution function of 
the number of contacted stations.at each step and the terminal relia- 


bility of the network. 
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111. COMPARATIVE ANALYSIS 


For the purpose of comparison, the computer program has been used 
to investigate a random communication network having 21 stations 
whose probability of contact is assumed constant. Three cases have 
been considered where the probability is equal respectively to 0.05, 
0.10, and 0.15. The sampling size is calculated to have an error of 


0.01 and the results are compared with the theoretical values. 


A. THE CALCUATION OF THE SAMPLING SIZE 
Suppose that one station tries to send the message to all other 


Stations. The expected number of contacted stations is given by 


о ое p (3.1) 


The variance is given by 


Geman = 1) pq 
where n = the size of the network 
p = the probability of contact 
q=1-p. 


Equation (3,1) becomes 


| 


in the case where n = 21 and p = 0.05. 
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Let Z(1) designate the random number of contacted stations at the 
first step in one sample. Since the population is infinite, the mean 
and the standard deviation of the sampling distribution of means are 


given by 


|2 у= 241) à - 


аа 972(1) ---- 
UN 
where N = the sampling size. 

For large values ОЁ М (М > 30) the sampling distribution of means 
is approximately a normal distribution with mean ы and standard de- 
viation oN . The accuracy of the approximation improves as N 
gets larger (Central limit theorem). 

Suppose that one would like the value of Z(1) differ from the 
population mean | by 0.01 with a confidence level of 95.45%, one has 


to satisfy the inequality 


| 21) -p | < 2071) < 0.01 (252) 
от „ №2 (п-1)ра < 10-2 
от м у 4.10% . (n-1)pq = 4.104 , (21-1) (0.05) (0.95) 


then NY 3.8 x 104 


The program has been run through 3.8 x 104 iterations. The results 
of the simulation are compared with the theoretical values in the fol- 


lowing section. 
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В. COMPARISON BETWEEN SAMPLING AND THEORETICAL VALUES 


a. Average Numbers of Contacted Stations at the First and Second 
Steps 


The average numbers of contacted stations are given by [5] 


Il 


Z(1) = (n-1)p 


Z(2) = (n Def: - (а-р2у 2, 


ТП 


They are calculated for three values of the probability p and compared 
With the results of the sampling in tables 3.1 and 3.2. The difference 


between sampling and theoretical results is small. 


Table 3:3. Values of Z(1). 


— ne 










Difference 





Probability Theoretical 





Sampling 


— me er ee ee oe 


0.9946 






| 
| 
| 
| 
| 
| 
| 
рез. 





Table 3.2. Values of Z(2). 


Se per АС ===> EEE о с 


| — = | 
Probability Theoretical Sampling Difference 








0.05 | 0.8825 0.8708 0.0117 
0.10 | 8912890 3.1334 0.0045 





43 





The sampling size has been calculated for the value of probability 
equal to 0.05. Table 3.1 shows that the difference between the sampling 
апа theoretical values 2(1)- 1 is 0.0022, less than 0.01 as it is 
expected. 
b. Average Total Numbers of Contacted Stations 
Fig. 3.1 gives the numbers of stations which will be eventually 
contacted for different values of probability. They are tabulated in 


Table 3.3 along with the corresponding theoretical values. 


Table 3.3. Average total number of contacted stations 


| 
Probability Theoretical | Sampling Difference 


——————— сызы 





TM 
-0.07 | 


The theoretical [5] and sampling curves giving the average numbers of 
contacted stations up to step K are similar. Both show that the 
propagation of message is over between steps 4 and 6. 
C. Terminal Reliability 
The terminal reliability curves in [5] and in Fig. 3.2 are 
very similar. They show that the highest reliability occurs between 
steps 2 and 3 and is negligible from step 6. The terminal reliabil- 


ity is the measure of the speed of propagation of the message. 
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Table 3.4. Values of highest terminal reliability 





a o 0o 


| | | 
Probability | Theoretical | Sampling ; Difference 


шн а — — — E 
0.10 | 0.175 | 0.175 | 0.0 
| ERUNT - 
0.15 | 0.34 ‚ 0.32 | 0.02 | 
NEU wc WI o 


d. Probability Distribution of the Number of New Stations 
Contacted at Each Step 


Fig. 5.3, 3.4, and 3.5 give the probability distributweonsert 
new stations contacted at each step, they correspond to the formula 
(1.3), they are very similar to the theoretical curves in [5]. 

In summary, the computer program has yielded results almost 
identical to the theoretical values if the sampling size is large, 
between 4 x 10% and 10%. Since the program was written without any 
constraint, any assumption on the probability of communication between 
pairs of stations, it could be used with confidence to study any ran- 
dom network on the condition that the matrix of probabilities is 
initialized to values which reflect each particular circumstance. As 
an example, in the next chapter, the program is used to investigate 


a communication network whose destruction is supposed normally dis- 


tributed around a nuclear explosion or an earthquake. 
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Fig. 3.1. Average message propagation 
x scale = order of step (2 units per inch) 
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Fig. 3.3. Probability distribution of the number of contacted 
stations at each step. 
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Fig. 3.4. Probability distribution of the number of contacted 
stations at each step. 


Probability = 0.1 
x scale = number of contacted stations (3 units per inch) 
y scale = probability distribution (0.1 unit per inch) 


k order of step 
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Fig. 3.5. Probability distribution of the numbers of contacted 
stations at each step. 


Probability = 0.15 


x scale = number of contacted stations (3 units per inch) 


y scale = probability distribution (0.1 unit per inch) 


k = order of step 





IV. INVESTIGATION OF A GENERAL MODEL 


А. REPRESENTATION OF THE MODEL 
Consider a communication network of n stations which may be rep- 


resented by an n x n matrix whose elements are 1 or 0 as 


1 273 n 
110 1 0 0 1 1 
210 O 1.0 0 0 
311 0 0 0 1 
0 
(4.1) 
Р- 1 
po 


Suppose there is a perturbation near station l, whose effect is to 
destroy the links of communication according to a normal distribution 
N(0,g-) around a center. This means that the destruction is maximum 
at the point of impact and decreases radially according to the 


equation 





x 2/_2 
р(х) ] 1 1/287 dt. 
O og Vår 
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Fig. 4.1. Illustrating the equation (4.2) 


where p(x) = probability of survival of a link of a station 
x = distance of a station to the impact 
o = standard deviation of the destruction 


The probability of communication between stations i and j is 


Pij = РО) p(xj) 
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Suppose that the stations are numbered according to increasing dis- 
tance and the network is complete, then all elements of the matrix P 
have to be replaced by the values of equation (4.2). The value of 
the standard deviation measures the extent of the destruction. The 
simulation was made for three cases where station 21 is one, two and 


three standard deviations away from the impact. 


first case g = 21 extensive destruction 
21 ' : 

second case oO” = q 10.5 medium destruction 

third case g-2- 7 minor destruction 


The subscripts of the stations which represent their distance to the 

impact, have been expressed in standard units and the subroutine NDTR 
has been used to initialize the matrix P. The sampling size has been 
taken equal to 3.8 x 104. The results are summarized in the following 


section. 


В. RESULTS OF THE SIMULATION 

When the farthest station is one standard deviation from the im- 
pact, i.e, the case when the destruction is extensive, the message is 
not carried to many stations. On the average, 1.26 stations receive 
the message, the link probabilities being very low. In the case of 
medium destruction, the farthest station two standard deviations from 


the impact, 4.17 stations receive the message. The propagation 
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ceases at step 5. 7.45 stations receive the message when station 21 
is three standard deviations away from the impact. The above results 
are obtained with the assumption the communication network forms a 
complete graph, 1.e., there is a link between any pair of stations. 
It is usually not the case in real life, therefore the above results 
are only the upper bound of number of stations which will receive the 
message in case of a war. 

If one compares with the model where the link probability is con- 
stant, one finds the latter model always gives much higher results 
than the current model which is nearer to reality. Table 4.1 summar- 


izes the results of two models. 


Table 4.1. Number of contacted stations in two models 
compared. 


DAA AAA 





A A e XA 





Probability — Mattel Standard Deviation . Normal Model 
0.05 | 4.40 21 АЧА. 
| | 
0.10 | 14.07 10.5 14.17 
! 


0.15 | Top) 7 | 7.45 


Fig. 4.2 shows the pe message propagation in the three cases. 

The weak connectivity of the normal model is also smaller than the 
Mattei model, it is equal to 0.06, 0.20 and 0.35 when the farthest 
station is away from the impact one, two and three standard deviations 
respectively. The link densities are respectively 0.08, 0.36 and 
0.67. All results point out that in the normal model, the propagation 
stops rapidly. The terminal reliability in Fig. 4.3 is a monotoni- 


cally decreasing function in the first case, attains the maximum 
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value at step 3 for the second and third cases. It means that the 
propagation increases, attains the culmination then decreases and 
stops in average at step 5. The curves of probability distribution 
function in Figs. 4.4, 4.5 and 4.6 of the number of contacted stations 


at each step show also small values compared with the Mattei model. 
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Fig. 4.2. 
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Average message propagation 
order of step (2 units per inch) 


number of contacted stations (3 units per inch) 
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Eig. 4.3. 


x scale 


y scale 





Terminal reliability 
order of step (1 unit per inch) 


terminal reliability (0.1 unit per inch) 
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Fig. 4.4. Probability distribution of the number of contacted 
stations at Gach step 


Standard deviation = 21 


x scale = number of contacted stations (3 units per inch) 


y scale = probability distribution (0.2 unit per inch) 
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order of step 
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Probability distribution of the number of contacted 


stations at each step 


Standard deviation = 10.5 


x scale 
y scale 


k 


number of contacted stations (3 units per inch) 


probability distribution (0.2 unit per inch) 


order of step 
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Fig. 4.6. Probability distribution of the number of contacted 
stations at each step 


Standard deviation = 7 


x scale = number of contacted stations (3 units per inch) 


y scale = probability distribution (0.2 unit per inch) 


k 


order of step 
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V. CONCLUSION 


The analysis of the message propagation through a random communi- 
cation network has been studied by various investigators. To be able 
to solve the problem, each has adopted different assumptions on the 
probability of communication between stations. The main effort of 
this paper has been to solve the same problem without having to make 
any simplifying assumption on the link probability. The latter could 
vary to fit a particular circumstance. 

A flowchart has been drawn and a program for its implementation 
on the digital computer has been written. The program gives the de- 
sired quantities of a random communication network whose link 
probabilities could be varied to fit a particular situation by the 
Monte Carlo technique. The program has been tested with the model 
whose link probabilities were assumed constant and whose characteris- 
tics could be derived analytically. The results were found in 
agreement with the analytical results. The program has also been 
used on a model where the destruction was supposed to be normally 
distributed around the point of attack. The results also agree with 
the intuitive idea that, as the effectiveness of the destruction of 
the network increases, the propagation of message stops rapidly. 

In the computer program, the time taken by the message to travel 
between any pair of stations is assumed constant, which is not true 
in real life. Some stations take longer time to process and send the 
information than the others. It is proposed for future work that the 
time tij; time taken by station i to process information and send it 


to station j, could be varied at will. With this second 
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generalization added, the simulation will give the average time of the 
propagation of the message and what routing procedure will give the 


fastest communication between various stations. 
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